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We discuss some properties of the generalized entropies, called Renyi entropies and their appli- 
cation to the case of continuous distributions. In particular it is shown that these measures of 
complexity can be divergent, however, their differences are free from these divergences thus en- 
abling them to be good candidates for the description of the extension and the shape of continuous 
distributions. We apply this formalism to the projection of wave functions onto the coherent state 
basis, i.e. to the Husimi representation. We also show how the localization properties of the Husimi 
distribution on average can be reconstructed from its marginal distributions that are calculated in 
position and momentum space in the case when the phase space has no structure, i.e. no classical 
limit can be defined. Numerical simulations on a one dimensional disordered system corroborate 
our expectations. 
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I. INTRODUCTION 

The measure of the extension of phase space distribu- 
tion of quantum states tells us important information on 
the degree of ergodicity and at the same time the degree 
of localization. These information are directly connected 
to the caoticity of the underlying classical dynamics if 
the latter is meaningful. In the ergodic regime trajecto- 
ries visit every corner of phase space hence the quantum 
states associated to such orbits are expected to be ex- 
tended. On the other hand regular islands trap classical 
trajectories and the corresponding states are localized. 
Therefore the extension properties of the eigenstates di- 
rectly reflect the nature of classical dynamics. For this 
purpose the Shannon entropy or information content has 
been used widely as a measure of complexity of quan- 
tum states. This and other generalized entropic mea- 
sures have been invoked by Zyczkowsky in Ref . [lj as mea- 
sures of chaotic signatures. In subsequent work 2 this idea 
has been elaborated further and a direct correspondence 
between the complexity of quantum states and the un- 
derlying dynamics has been demonstrated in particular 
by projecting the quantum states onto a coherent state 
basis. Recently other works have showed the renewed 
interest in this fie\d£££&. We have to emphasize, how- 
ever, that even systems without classical limit, e.g. dis- 
ordered systems, have been involved in such phase space 
studies 7 - 8-9 . This latter topic is the main motivation of 
our present work as well. 

In the present work we give further arguments in favor 
of the application of generalized entropies, the Renyi en- 
tropies, for the characterization of quantum phase space 
distributions, however, we point out some problems in 
connection to the calculation of these parameters for con- 
tinuous distributions. As a remedy for these problems we 
show that the differences of Renyi entropies are on the 
other hand free from these anomalies. 



Furthermore we will show that indeed these entropic 
measures give important information concerning the lo- 
calization properties of phase space distributions espe- 
cially the Husimi distribution. We will also give argu- 
ments and numerical proofs that in fact in the case of 
wave functions of disordered systems there is no need to 
calculate the Husimi distributions themselves especially 
because the average localization properties of the Husimi 
functions of a set of states can be obtained from the aver- 
age properties of the marginal distributions of the Husimi 
functions. 

Using different techniques similar results have been ob- 
tained for a particular quantity, the participation ratio in 
Refs. 00. Our approach, however, is more general. 

In the next section we introduce the basic ideas and 
tools that have been widely used in the literature, namely 
the participation number (ratio) and the Shannon en- 
tropy for the case of discrete distributions. We also 
show that these quantities give roughly the same infor- 
mation, however, using these parameters a new quantity, 
the structural entropy, can be defined that contains im- 
portant information concerning the shape of a distribu- 
tion. It is also shown that these parameters are nothing- 
else but some special cases of the differences of Renyi en- 
tropies. Section II contains merely the revision of what 
has been published before and we conclude this section 
analyzing the problem of continuous distributions and 
showing that the above mentioned differences of Renyi 
entropies are free from the divergences. In Section III 
we elaborate the appropriate generalization of these pa- 
rameters for continuous distributions. In Section IV we 
introduce the Husimi representation of quantum states 
and describe some of its properties. In Section V Renyi 
entropies are applied for Husimi distributions and it is 
shown that the properties of its marginal distributions al- 
ready give a qualitative picture that for special cases may 
quantitatively be correct, as well. In Sect. VI numeri- 
cal simulations for the one dimensional Anderson model 
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provide important verification of the results presented in 
Sect. V. Finally some concluding remark are left for Sect. 
VII. 



II. BASIC IDEAS 



which has been termed as the participation ratio in the 
literature. These two quantities satisfy the following 
inequalities^ 

< q < 1 (7a) 
0< S str <-lnq. (7b) 



The extension of a discrete distribution of a state is 
often characterized by its entropy, S, or by its second 
moment, D. Both of these quantities, i.e. exp(S) and D 
practically measure the same thing, namely the number 
of amplitudes that mainly contribute to the expansion of 
the state over a suitable basis. 

Let us consider a wave function 'J that is represented 
by its expansion over a complete basis set on a finite grid 
of N states <fii 

N N 

$ = J2a<t>i with J]|c 4 | 2 = i (l) 

»=i i=i 

Note that each of the coefficients, Qi = \ci\ 2 , obey the 
condition 



< Qi < 1, 



(2) 



and they sum up to unity. Then the usual definitions of 
participation number, D, and entropy, S are 



JY 



and 



S = - Qj In Qj 



(3) 



The parameter D tells us how many of the numbers Qi 
are significantly larger than zero. For example if only one 
of them is unity and the rest is zero, then D = 1. Oth- 
erwise if Qi = 1/N homogeneously, then we get D = N. 
Similar properties can be shown to hold for exp S 1 , there- 
fore it is easy to show that the two quantities provide 
roughly the same information 



In A 



(4) 



i.e. both, S and lnD describe the extension of the dis- 
crete distribution. That is the reason for calling D as 
the number of principal components. The close relation 
between S and D, Eq. QJ, has often been overlooked and 



presente' 



10.11 



as an interesting similarity. However, as it 



has been demonstrated in Ref. [12j and applied in several 
studies later 12 i 13 i 14 i 15 i 16 , the difference 



S str = S -InD 



(5) 



is a meaningful and most importantly a nonnegative 
quantity that turns out to be very useful in the char- 
acterization of the shape of the distribution of the prob- 
abilities Qi. That is the reason why it has been termed as 
structural entropy of a distribution. Moreover the value 
of the participation number normalized to the number of 
available components is an important partner quantity 



D 

TV' 



(G) 



Generalized entropies have been introduced by Renyiii 
in the form of 



Rm — 



1 N 

1 — m ^-^ 



(8) 



that monotonously decreases for increasing m. For the 
special cases of m = 0, 1, and 2 we recover the total num- 
ber of components, the Shannon-entropy and the partic- 
ipation number 

R = IniV lim R m = S, R 2 = InD. (9) 

m — >1 

Notice that the order m in JSJ is not necessarily integer. 
We can readily realize that the parameters, JSJ and iJBJ, 
are nothing else but the differences^ 3 - of the special cases 
of i?o, R\, and R 2 , i.e. 

S str = Ri-Ri (10a) 
-bag = R -R 2 - (10b) 

A number of applications have been presented to datei 3 - 
showing their diverse applicability starting from quan- 
tum chemistry— up to localization in disordered and 
quasiperiodic systems 1 ^ up to the statistical analysis of 
spectrai 6 .. 

In the present work we are going to extend this for- 
malism to continuous distributions and show again that 
the differences of Renyi entropies are good candidates for 
the characterization of them. 

The problem with a continuous distribution p(x) is 
that even though normalization requires 



dxp(x) = 1, 



(11) 



it is clear that p{x) is a density of probabilities, therefore 
it is the quantity p(x) Ax, the probability associated with 
the interval [x, x + Ax] that is restricted to the [0, 1] in- 
terval and not the value of p(x) itself. Hence even though 
we may always expect p(x) > the condition p(x) < 1 is 
generally not fulfilled. 

Then the obvious generalization of the Renyi entropies 
JHl for normalized continuous distributions would read 



as 



Rr, 



1 



In / dx [p(x)] r 



(12) 



Hence the definitions of the participation number and the 
entropy © would look as 



D~ 



S = 



dx [p(x)Y 



dxp(x) hip(x). 



(13) 
(14) 
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Let us apply this definition to a Gaussian distribution 
with zero mean and variance a 



p(x) 



1 



V2 



: exp 



x 



(15) 



The formulas derived using this p(x) will be useful when 
applying for the problem of Husimi distributions later 
due to the Gaussian smearing contained in those phase 
space functions. Putting Eq. I|15|) in Eq. (|12f> we obtain 
for m > 



hi 



m — 1 



ln(V27rer), 



(16) 



which tells us that as a — > or a — > oo the Renyi en- 
tropies diverge. However, they do that uniformly, i.e. 
independently of m, hence their differences remain finite. 
This is an important advantage of our formulation that 
will be elaborated further in the subsequent part. In par- 
ticular in the case of the Gaussian (|15fl we obtain q = 
and 



i?2 



1 



-(1 - In 2) = 0.1534. 



(17) 



This is the value of the structural entropy that describes 
a Gaussian distribution in one dimension. 

Next let us consider the above Gaussian on a finite 
interval, —L/2 < x < L/2 and assume that beyond this 
interval p(x) = 0. This construction allows us to study 
how the limit of l|l(j|) or (|17|) is approached as for fixed a 
the interval tends to infinity (or for a fixed L the width 
a — ► 0). To this end the normalization is taken over the 
finite interval [—L/2, L/2] so the distribution function 
(II 51) should be modified as 



p{x) 



1 



-x 2 /2a 2 



(18) 



where $(cc) denotes the error function and the scaling 
parameter £ = L/a has been introduced. The Renyi 
entropies will depend on L and a as 



R m (L,a) = R m ((X),a) 



1 



1 



In 



(19) 



where -R m (co,er) is given in Eq. 116(1 . Again the m— 
independent — In a divergence appears. Turning to the 
special cases of q and S s t r as deduced from i?o, Ri, 
and i?2 using Eqs. flT)j) as before we find that they are 
uniquely determined by £ and are free from this type of 
divergence. In particular since Rq = InL, 



?(0 



20F mi V8)f 
£ *(£/2) 



and 



S s tr(0 — Sg 



£e 



^7T^((/V8) 



In 



$(£/2) 



(20) 



, (21) 



where Sf tr is given in Eq. I|17(l . The participation ratio, 
q(^), in the limit £ — > oo (i — > oo for fixed it or er — > for 
fixed L) tends to zero as q(£) w 2^/tt/^ while S str (0 — > 
S'fj,.. On the other hand in the other limit of £ — > (<r — > 
oo for fixed L or L — > for fixed a) we see that g(£) w 
(l-£ 4 /720) and SWCO ~ £ 4 /1440 therefore the relation 
S str w (1 - g)/2 is also fulfilled^. We would like to 
emphasize that no divergences are found for parameters 
q and S'str and they show well defined behavior in either 
limit. 



III. COARSE GRAINING 

Now let us turn to a more general investigation of our 
parameters. In this section we provide a natural gener- 
alization of the calculation of the parameters q and S s t r 
for continuous distributions. 

Let us consider a disjoint division of the interval, O, 
over which the distribution p(x) is defined. Each of these 
subintcrvals have an index i running from 1 to N and a 
size of u>i, such that = ^- Then let us define a 

characteristic function Xi( x ) 



X G LOi 

otherwise 



These functions are orthogonal, 

dxXi{x)Xj(x) = uJiSij, 



(22) 



(23) 



where &y is the Kronecker delta. The coarse grained 
value of the distribution p(x) in interval i is p(x)u>i = Qi, 
if x e u>i, more precisely 



(24) 



Qi= I dxp(x)xi{x). 



This way our coarse grained approximation to the density 
function is 



p{ X ) = Y.-^ ( 25 ) 

I 

that obviously satisfies the normalization condition 

/ dxp(x) = ^Q t = l. (26) 

i 

On the other hand the integral of the square of this func- 
tion using Eqs. lj23*|) and is 



J t-r 1 uji uj s— 1 



(27) 



For sake of simplicity we have used (and will use from 
now on) an equipartition, a;, = uj = Q/N. Now let us 
calculate the participation number, D 7 and entropy, S, 
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using Eqs. I|13fl from the discrete sums over the probabil- 
ities Qi. First it is clear from Eq. (|27[l that 

- In D = In IpO)] 2 + In (28) 

Similar procedure leads to 

S = — ^^QilnQi = — J dxp(x) lnp(x) — Inu, (29) 

i 

The appearance of the term In uj in both of these equa- 
tions shows another type of divergence originating from 
the subdivision of the interval f2, since the limit of 
N — > oo corresponds to lu — > 0. Therefore a naive ap- 
plication of these quantities may encounter severe con- 
ceptual and also numerical difficulties depending on the 
value of uj. On the other hand we may conclude once 
again that the parameters q and S s t r are free from this 
type of divergence, as well as using Eqs. JSJ, ||SJ, iJBJ 



ft J dx [p{a 



(30a) 



dxp(x) hxp(x) + In / dx [p(x)] 2 {30b) 



In q = In 



This way we have shown how to apply the formalism 
developed for discrete sums for the problem of continuous 
distributions. 



IV. PHASE SPACE REPRESENTATION OF 
QUANTUM STATES 

One of the most well-known phase space distributions 
that is widely applied in statistical physics is the Wigner- 
function associated with the quantum stated tpi x ) 



W(x,p) 



dx' e~ lpx ' /h ip* [x 



if) [ x + 



( 31 ) 

From now on p denotes momentum and for sake of sim- 
plicity we consider only one degree of freedom resulting 
in a two dimensional phase space of [x,p). 

It is known that W(x,p) is bilinear, real and for a 
complete orthonormal set of ipi functions the correspond- 
ing Wigner transforms also form a complete orthonormal 
set— The marginal distributions of W(x,p) have an im- 
portant physical meaning 



dpW[x,p) = |V(x)| 2 , 
dxW(x,p) = |0(p)| 2 , 

where (pip) denotes the Fourier transform of ipi x ) 
1 



<Pip) 



dxiP{x)e- ipx/h . 



(32a) 
(32b) 

(33) 



The only major disadvantage of the Wigner distribution 
is that is may attain negative values albeit in a region of 
phase space smaller than h. It has been shown already 
by Wigner— that there exists no phase space distribution 
that would have all the above properties and besides that 
to be nonnegative. 

Another very popular phase space distribution is the 
Husimi function^ that is obtained as the Gaussian 
smearing of the Wigner function, W(x,p) 



p(x,p) 



exp 



dx'dp' W(x',p') x 

ix-x') 2 (p-p 1 ) 2 



2(7? 



(34) 



where a x a p = H/2 ensures minimum uncertainty. 

It is knowni& that the Husimi function in bilinear, real 
valued and nonnegative but unfortunately produces an 
over-complete set of functions and moreover the marginal 
distributions do not have such a transparent meaning as 
in Eq. I|32|) . In fact the latter point can be refined. Let 
us calculate these marginal distributions and will find 
that indeed they are the Gaussian smeared distributions 
in x— and p— representations, respectively^. In order to 
show this we write (|34ll in the form of a convolution of 
W(x,p) with two Gaussian functions, g aw (x) and g ap {p) 
of the form of Eq. (|15|) 



p(x,p) = J dx'dp' g ax {x-x')g ap (p-p')W(x',p'). (35) 
Then similarly to Eq. 



dpp(x,p) = CO), 
dxp(x,p) = i](p), 



(36a) 
(36b) 



where the marginal distributions are nothing else but 
smeared distribution obtained from quantum state, tp( x ) 
and (p{p), respectively 



Cix) = I dx' 9( r x (x-x')\i>ix')\ 2 , (37a) 
V(p) = / dp'g^ip-p'Mip')] 2 , (37b) 



It is clear from the definition of the Husimi distribution 
that it is normalized as 



Jdxd PP ix, P ) = l, 



(38) 



therefore the smeared a:— representation of ip(x), Ci x ) an d 
the smeared p-representation of (pip), r](p) are normalized 



as 



J dx(ix) = J dpriip) = 1. 



(39) 
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To complete this section we mention that the Husimi 
representation of a quantum state ip(x) is nothing else 
but its projection onto (i.e. the overlap with) a coherent 
state with minimal uncertaint y 4 ' 5 ! 8 ! 21 j3 X:P (x') 



pj,(x,p) = \(P(x,p)\i/j)\ 2 = 



dx'p* XtP (x')^(x') 



(40) 



the coherent state is a Gaussian centered around the 
phase space point (x,p) 



1 



2n<r 



1/4 

j i ex P 



{x' - xf 

4<7 2 



ipx' jh 



case is the Husimi distribution, see Eq. (|34() or Eq. i|40|l . 
Its information content, or Shannon entropy 22 <|43[) can 
be related to the Shannon entropy of the marginal distri- 
butions, S[Q and S[rj\ defined in l|14|) for ((x) and r)(p) 
(|37|l . respectively. Let us note that the Husimi distribu- 
tion can be written in the form 



where 



p(x,P) = ((x)ij(p) + 5(x,p), 



dxS(x,p) — / dpS(x,p) = 0. 



(46) 



(47) 



(41) The Shannon entropy then obbeys^ the following relation 



V. RENYI ENTROPIES OF PHASE SPACE 
DISTRIBUTIONS 

In this section we describe how to characterize localiza- 
tion or ergodicity using the ingredients explained in the 
previous sections: (i) Husimi representation of the quan- 
tum states ijj and (ii) the Renyi entropies, especially their 
differences. 

First of all let us introduce the Renyi entropies of the 
Husimi function. In analogy with the definition l|12|) 



S[p]+\nh = S[(] + S[rj\ + SS 



(48) 



1 



1 — m 



1, / ^[hp(x lP )] m (42) 



which now contains the arbitrary parameter h, that nat- 
urally behaves as the minimum possible volume provided 
by the Heisenberg uncertainty principle, i.e. it should be 
chosen as the Planck's constant. We have to note that 
each R m contains In h additively, that diverges in the 
classical limit h — *■ but drops out when differences of 
the entropies are taken. Definition (|42l) for a compact 
phase space of volume f2, for instance, provides for the 
special case, Rq — ln(Q/h), i.e. it measures the size of 
the full phase space in units of h. Furthermore 

Ri = S = — J dxdp p(x,p)\n[hp(x,p)] (43) 

i? 2 = -ln^y dxdp[p(x,p)] 2 ^j . (44) 

These are in accordance with Boltzmann's original def- 
inition, since for a distribution that is constant over a 
volume r < Q and zero otherwise we obtain 



S = ln(r//i), 



and 



q = r/n, 



(45) 



i.e. S measures the size of phase space, T, where p is 
nonzero in units of h and q measures the portion of phase 
space where p is different from zero. 

In order to relate the entropy of the total Husimi func- 
tion to that of the marginal distributions let us invoke an 
important relation that has been proven for the Shannon 
entropy, S. Consider a distribution p(x,p) which in our 



where SS < 0. Equality is achieved if S(x, k) = ev- 
erywhere. This statement can be generalized to Renyi 
entropies where 



R m [p] + In h = R m [C] + R m M + 5Rr, 



(49) 



with 8R m < 0. Unfortunately there is no general law 
for the size of the 5R m , however, for the differences of 
Renyi entropies it may become only a small correction 
if 5(x,p) <ti C(x)ri(p). Furthermore for the parameters 
— In q = R 2 — Rq and S str = R% — R2 we should have cor- 
rections of 8R2 — SRo and 5R\ ~ <5i?2- These differences 
especially after averaging over several wave functions can 
be neglected, thereofore we arrive at the following ap- 
proximate relation for the average values of — In q and 

Sstr 

-\nq[p] » -lng[C]-ln«fo] (50a) 

S s tr[p} ~ S str [(] + SstrM (50b) 

Such relations reduce the calculations considerably as 
there would be no need to calculate the Husimi func- 
tions themselves and then calculate the Renyi entropies 
thereof. In a numerical application we will show below 
that indeed these relations do hold with small error. 

We have to stress that these approximate relations may 
hold apart from the trivial case of the wave packet (pre- 
sented next) for the average properties of a suitably cho- 
sen set of states and most importantly in the case of the 
lack of underlying classical dynamics. These limitations 
reduce its applicability for the investigation of the states 
of disordered systems, which is nevertheless the main aim 
of our study. 

As a simple example we elaborate the case of a distri- 
bution that is a Gaussian in both coordinates x and p 



p(x,p) 



1 exp (-x 2 /2a 2 - 2(ap/hf) 



(51) 



$(a/V8)$(\/27r/?) 

and normalized over the phase space bounded as —L/2 < 
x < L/2 and —nh/a < p < irh/a, where two cutoff length 



scales, L and a, have been introduced. Therefore the 
volume of the phase space will be Q — hL/a = 7/1. The 
ratios of the cutoff scales to the spreading width a yield 
the two dimcnsionless parameters in (|51|) . a — L/a and 
/3 = a /a. In terms of these quantities we obtain the 
relation 7 = a/3, which counts the number of cells of size 
h in phase space. 

Since in (|51|1 a x = a and a v = h/2a, the uncertainty 
relation, a x a v = h/2, is fulfilled. This p(x,p) is in fact 
the Husimi distribution of a real space Gaussian wave 
packet and it is a product of the limit distributions as 
given in 146|) with S(x,p) — 0. Consequently the q and 
S s tr values of the corresponding limit distributions, Ci x ) 
and T](p), obey the additivity property (|5U|l exactly. Dis- 
tribution C( x )j f° r instance, is obtained using (|36aj) and 
yields Eq. JT&J. Its q[(] and S s tr[C] parameters are given 
in Eqs. (|20|l and (|21l) . A straightforward calculation 
yields q[rj\ and S s tr[rj\, as well. 

Putting the Husimi distribution (frTTJl into we find 



0.4 



R m (a,(3) 



lnm 
m — 1 
1 

m - 



-ln2 



1 



\ [$(a/V8W0nV2)) n f 



This expression correctly yields R a — ln(a/3) = ln(7), the 
log of the number of cells of size h. Through relations 
HI U| l parameters q and S s tr can be obtained. Equivalently 
using the definition (|5I|I and (|30(l generalized for the case 
of the Husimi distribution like in 1|42JI and l|43|l we obtain 
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FIG. 1: The parameters S 3 tr (a, semi-log plot) and q (b, 
log-log plot) for a Gaussian wave packet as a function of the 
parameter t = a/oto, where a = L/a, L is the system size 
and a is the spreading of the Gaussian. Qo = s/Any. The 
different curves are parametrized according to 7, the number 
of cells within h. 



variable t — a/ao we obtain 



q(a,P) 
S str (a,P) 



Ma/V8Wf3V2ir)} 
a/3$(a/2)$(27r/3) 
I - ln(2) 



(53a) 



ae 



-a 2 /8 



In 



<S>(V2tt{3) V8n$(a/V8) 
$(27r/3)$(a/2) 



.$(v^7r/3)$(a/V8) 



(53b) 



There are a number of remarkable features of Eqs. I|52|) 
and i|53b[) . None of them contain Planck's constant, h, 
explicitly, however, the scaling variables a and /3 can- 
not be chosen independently as their product is just the 
number of cells of size h. Therefore let us keep a as a 
running variable and parametrize the functions with 7. 
Furthermore, let us note that by keeping 7 fixed the q(a) 
and S str (en) functions are symmetrical about a — ^/Any 
on a logarithmic scale. Therefore rewriting l|53|) in the 



,(7) 



(*) 



1 [$(d/t)$(cit)] a 
7 $(c 2 /i)$(c 2 i) 

1 - ln(2) 

/7 ,' te -( C it) 2 e -(ci/t) 2 ' 

2 j *( Cl t) + t$( Cl A) 

$(c 2 /t)$(c 2 t) ' 




hi 



$(ci/i)*(cii) 



(54) 



where ci = y / 7T7/ 7 2 and c 2 = \/2ci. These functions are 
plotted in Figure^ At t — 1 function q(t) is maximal and 
its value decreases with 7 as q(l) cx 7 _1 . Only the physi- 
cally relevant, 7 > 1 are plotted. The participation ratio 
has some very nice, simple behavior, q^(t) — > (^y^) -1 
for i — > co and q^(t) — > t/^Py for i — > 0. On the other 
hand in the same limits 5 s t r — » independently from 7 
showing that these limits correspond to one dimensional 
Gaussians in x(p-)directions for t — > oo(t — > 0). It can 
also be viewed as if a squeezing parameter made the dis- 
tribution more coordinate-like (momentum-like For 



intermediate values of t, i.e. if 1/^/7 < t < y/^y with 
7>1 the participation ratio is, q w 7 _1 and 5 s t r » 2iS^ r 
indicating that this distribution is a Gaussian in both di- 
mensions, x and p. This is the regime where the Husimi- 
function is a Gaussian in both directions and therefore 
shows a two-dimensional character. Both curves, q anf 
S s tr are symmetrical about t = 1 on a logarithmic scale 
of t, which is a direct consequence of the geometry of the 
phase space. 



VI. APPLICATION TO DISORDERED 
SYSTEMS 



Now we calculate Husimi functions of the eigenstates of 
a disordered one dimensional system. To be more precise 
we use a tight binding modelSi 

H = £n\n){n\ + (!»)<« + 1| + |n + l)(n|) (55) 



where V = 1 is set as the unit of energy and e„ are 
random numbers distributed uniformly over the interval 
[-W/2, . . . , W/2], where W characterizes the strength of 
disorder. Such a model has been investigated in phase 
space in Refs. □ and |3- The Husimi functions are cal- 
culated using Eqs. (f4*0")) and ijlT)) from the eigenstates of 
(|55p. The participation ratio, q, and the structural en- 
tropy, S s tr, are calculated according to Eq. i|30|) . We 
also calculated the Fourier transforms of the eigenstates 
and obtained smeared distributions according to Eq. Ij37(l 
both in real and Fourier space. Periodic boundary con- 
ditions were considered using L = 512. The phase space 
extends over — L/2 < x < L/2 and — 7r < p < it, its area 
is SI = 2-kL (H = 1 and the lower cut-off scale, the lattice 
spacing is set to unity, a = 1). Averaging is done over 
the middle half of the band. In fact, as pointed out by 
Ref. |a, as well, there is no need to average over many 
realizations of the disordered potential. 

When the full Husimi distributions of all states are cal- 
culated the computational time grows with L 4 . However, 
using the approximate relation of Eq. (|50J) . it reduces to 
roughly L 3 . This is obviously a considerable gain and is 
comparable to the one achieved in Refs. 00. 

The results are reported in Fig. |2J Here we have 
plotted the behavior of parameters l|5Uj l versus disorder 
strength, W , and compared to the approximate values 
obtained using (|50|l . The region where the most impor- 
tant variations of q and S s t r take place is W\ < W < 
W2, where the localization len gth ma tches the systems 
sizQ 24 J 25 ) x « L, i.e. Wi w ^JWb/L = 0.453 or' its in- 
verse approaches the size of phase space in p-direction, 
2ttA w 1, i.e. W 2 « V2tt10E = 25.68 

Let us analyze the expectations for the limiting cases 
of W — > and W — * 00. It is easy to see that the 
eigenstates for vanishing disorder are plane waves whose 
Fourier transform is a Dirac-delta (in fact two, due to the 
symmetry of the — p and p states). In a 'smeared' rep- 
resentation we obtain two Gaussians in p-representation. 




j ' 1 ' 1 ' 1 ' 1 



0.01 0.1 1 10 100 1000 
disorder (W/V) 

FIG. 2: Structural entropy, S s tr (a, semi-log plot) and par- 
ticipation ratio, q (b, log-log plot) as a function of disor- 
der, W in units of V, for a one dimensional Anderson model 
with L — 512. The squares stand for the states smeared in 
x-representation, the triangles for the states smeared in p- 
representation, the circles are calculated according to Eqs. 
I|5()|l . The dotted lines are simply guides for the eyes. The 
solid curve corresponds to the values of S 3 tr and q for the 
states in Husimi representation. The circles and the solid 
curve differ only a little. 



Using Eq. JIT} and (|30a|l we obtain that q = 2/ \fZ for 
the states both in p- and Husimi representation. Due to 
the Gaussian smearing in this limit the structural entropy 
attains its value of iS^ r as given in l(T7f). The other limit 
of W — > 00 is very similar. In that case the eigenstate 
in x-representation has a Dirac-delta character that is 
smeared to a Gaussian. This results in q = l/\/L and 
again a value of S s tr(W) — * Sf tr . All these limiting cases 
are recovered in Fig. [21 The figure shows that the ap- 
proximation l|50|l works very well. It is clear that the x-, 
p- and Husimi representations are therefore linked very 
simply. 



VII. CONCLUDING REMARKS 

In this paper we have presented some important re- 
sults concerning the applicability of Renyi entropies for 
the characterization of localization or ergodicity in phase 
space using the Husimi representation of the quantum 
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states. In fact it has been shown that the differences 
of Renyi entropies are free from those divergences that 
would naturally arise due to their application on contin- 
uous distributions. 

The marginal distributions of the Husimi function are 
pointed out to have important properties and simple con- 
nection to the states in x- and ^-representations. 

We have also shown numerically that for disordered 
systems the limiting distributions of the Husimi function 
provide most of the information that is needed to describe 
the Husimi functions themselves. Figure [21 provides a 
good demonstration of the duality between the x- and 
p-representations transparently. A detailed study over 



the Anderson model in one dimension and the Harper 
models are left for forthcoming publications. 
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